Preprint typeset in JINST style - HYPER VERSION 



Performances of a large mass ZnSe bolometer to 
search for rare events 



J.W. Beeman", F. Bellini^, L. Cardani* * N. Casa\\ de , I. Dafinei f , S. Di Domizio' 1 , 

F. Ferroni^, L. Gironi /g , A. Giuliani 1 , S. Nagorny e , F. Orio f , L. Pattavina e , 

G. Pessina^, G. Piperno foc , S. Pirro^, E. Previtali g C. Rusconi^, C. Tomei c , M. Vignati c 

"Lawrence Berkeley National Laboratory , Berkeley, California 94720, USA 
b Dipartimento di Fisica - Universitd di Roma La Sapienza, 1-00185 Roma - Italy 
C INFN - Sezione di Roma, 1-00185 Roma - Italy 

" Dipartimento di Scienze Fisiche e Chimiche - Universitd degli studi dell 'Aquila, 1-67100 
Coppito (AQ) - Italy 

e INFN - Laboratori Nazionali del Gran Sasso, 1-67010 Assergi (AQ) - Italy 
J Dipartimento di Fisica - Universitd di Milano Bicocca, 1-20126 Milano - Italy 
g INFN - Sezione di Milano Bicocca, 1-20126 Milano - Italy 
INFN - Sezione di Genova, 1-16146 Genova - Italy 

'CSNSM, Centre de Spectrometrie Nucleaire et de Spectrometrie de Masse, Bdtiment 108, 
Campus d'Orsay 91405 Orsay, France 



E-mail: laura . cardani@romal . inf n . it 



ABSTRACT: Scintillating bolometers of ZnSe are the baseline choice of the LUCIFER experiment, 
whose aim is to observe the neutrinoless double beta decay of 82 Se. The independent read-out of 
the heat and scintillation signals allows to identify and reject a particle interactions, the dominant 
background source for bolometric detectors. In this paper we report the performances of a ZnSe 
crystal operated within the LUCIFER R&D. We measured the scintillation yield, the energy res- 
olution and the background in the energy region where the signal from OvDBD decay of 82 Se is 
expected with an exposure of 9.4 kg -days. With a newly developed analysis algorithm we im- 
proved the rejection of a events, and we estimated the increase in energy resolution obtained by 
the combination of the heat and light signals. For the first time we measured the light emitted by 
nuclear recoils, and found it to be compatible with zero. We conclude that the discrimination of 
nuclear recoils from j3/y interactions in the WIMPs energy region is possible, but low-noise light 
detectors are needed. 

KEYWORDS: Bolometer, Neutrinoless double beta decay, Dark Matter. 



* Corresponding author. 



Contents 

|l| Introduction I 

|2[ Experimental setup [§ 

g Scintillation | 

f|, Energy Resolution |] 

|5[ Pulse shape discrimination § 

g Background 

[7[ Conclusions 



1. Introduction 

Bolometers are solid state detectors in which the energy release coming from particle interactions 
is converted to heat and measured via their rise in temperature. They can provide excellent energy 
resolution and low background, and are used in particle physics experiments searching for rare 
processes, such as neutrinoless double beta decay (OvDBD) and Dark Matter interactions. 

In the last years, an extensive R&D allowed to increase the mass of bolometric detectors and, 
at the same time, to understand and reduce the sources of background. The largest bolometric 
detector operated up to now, Cuoricino, was made of about 40.7 kg of Te02 and took data from 
2003 to 2008, demonstrating the potential of this technique [jlj The evolution of Cuoricino, 
CUORE, will search for the OvDBD of 130 Te [g, J§] using an array of 988 Te0 2 bolometers of 750 g 
each. Operated at a temperature of about 10 mK, these detectors provide an energy resolution of 
a few keV over their energy range, extending from a few keV up to several MeV. The measured 
resolution at the Q-value of the decay {Q = 2527 keV ^) is about 5 keVFWHM; together with 
the low background and the large mass of the experiment, this will provide a 1 a sensitivity to the 
OvDBD of 130 Te of the order of 10 26 years. 

To further increase the sensitivity, an intense R&D is being pursued to lower the background in 
the OvDBD region. The main source of background is due to a particles, coming from radioactive 
contaminations of the materials facing the bolometers [Q]. A way to discriminate this background 
is to use a scintillating bolometer [g]. In such a device the simultaneous and independent read- 
out of the heat and the scintillation light permits to discriminate events due to /3/y, neutron and 
a interactions thanks to their different scintillation properties. Unfortunately Te02 crystals do 
not scintillate and different compounds are being studied. Among these are CdW04 (OvDBD 
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candidate 116 Cd, Q = 2814 keV @]), ZnSe [g] ( 82 Se, Q = 2997 keV [|T(J]), and ZnMo0 4 [JTT], g 
DL ( 100 Mo, 2 = 3034 keV rtTl). 



In this paper we investigate the performances of a ZnSe scintillating crystal for the OvDBD 
search in terms of energy resolution, capability of discriminating a particles and internal radioac- 
tive contaminations. We also study the light emitted by nuclear recoils and j8 fy interactions at 
100 keV, discussing the possibility of using ZnSe bolometers as Dark Matter detectors. 

2. Experimental setup 

The data here presented come from a series of runs performed at Laboratori Nazionali del Gran 
Sasso (LNGS) in Italy, inside the CUORE R&D facility. The bolometer under study is a 431 g 
ZnSe crystal with cylindrical shape (height 44.3 mm and diameter 48.5 mm). The synthesis of 
the ZnSe powder, including purification and formatting for crystal growth [TQ 17], was made in 
ultra-clean fused quartz reactors at SmiLab Svitlovodsk (Ukraine). 

To detect the scintillation light, many different light detectors (LD) were faced to the ZnSe 
crystal. They consist of pure Germanium slabs of 50 mm diameter and with variable thickness 
(300 - 600 jUm) and they were operated as bolometers to obtain good performances at cryogenic 



temperatures Q18Q. We covered a face of the LD with a thin layer of Si02 (60 nm) and obtained an 
increase of the light absorption by ~ 16 - 20%, as already observed in Ref. [Q|. Moreover, we 
measured an increase of ~ 25% in the collection efficiency by surrounding the crystal with a 3M 
VM2002 reflecting foil. 

The temperature sensors of both the ZnSe and LD were Neutron Transmutation Doped (NTD) 



Germanium thermistors [20|, coupled to the ZnSe and Ge surfaces by means of epoxy glue spots. 
For redundancy, the ZnSe crystal was equipped with two sensors, in the following referred as 
ZnSe-L and ZnSe-R. The detectors were held in a copper structure by Teflon (PTFE) supports 
and thermally coupled to the mixing chamber of a dilution refrigerator which kept the system at a 
temperature around 10 mK. To read the signals, the thermistors were biased with a constant current. 
The voltage signals, amplified and filtered by means of an anti-aliasing 6-pole active Bessel filter 
(120 dB/decade), were fed into a NI PXI-6284 18-bits ADC operating at a sampling frequency of 
2 kHz. The Bessel cutoff was set at 120 Hz for the LDs, and at 70 Hz for the ZnSe. Further details 



on the cryogenic facility and the electronic read-out can be found in Refs. J2 1 , |22|, [23]]. 

The trigger was software generated on each bolometer. When it fired, waveforms 5 s long on 
the ZnSe and 250 ms long on the LD were digitized and saved on disk. In addition, for every trigger 
occurring on the ZnSe, the LD waveform was acquired irrespective of its trigger. The off-line 
analysis computes the pulse height as well as pulse shape parameters based on the optimum filter 



algorithm [ |24[ , [25| ]. In addition, the amplitude of the light signals is computed with an algorithm 
that allows to lower the energy threshold of the LD using the knowledge of the signal time delay 
with respect to the ZnSe (see details in Ref. [|26|]). 

Compared to Te02 crystals, we observed that the ZnSe bolometer cooled very slowly, likely 
because of a heat capacitance excess at low temperature. Indeed, despite of the fact that the cryostat 
temperature was kept constant, the temperature drift of the ZnSe was large enough to influence the 
detector response still after 30 days of data-taking (Fig. [I]). The resistance of the thermistor, which 
has a steep dependence on the temperature, increased of about 6% during this data acquisition 
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Figure 1. Change in the thermistor response induced by the slow cooling of the ZnSe bolometer, whose 
temperature decreased of about 1% during the first 700 hours of measurement. Left: increase of the thermis- 
tor resistance due to the cooling of the crystal. Right: amplitude of pulses with energy of about 1460 keV 
as a function of time. In the last 150 hours a calibration measurement was performed, resulting in a higher 
rate. 



period, resulting in a corresponding increase in the signal height. This behavior was also observed 
in bolometers made of BGO crystals [27]. 

The calibration of the ZnSe pulse amplitude was performed by means of a 228 Th y-source 
placed inside the cryostat external lead shield. We observed that the calibration function derived 
from 7 peaks was not applicable to energy deposits induced by a particles, which were shifted 
by +22% with respect to their nominal energy. For this reason, in the following we will refer to 
the energy estimated from the y calibration as "Energy ee " (ys and electrons interact in the same 
way in the bolometer). When the calibration is derived from a peaks produced by the internal 
contaminations, the energy units will be indicated as "Energy^". 

To evaluate the discrimination power between /j /y and a events in the energy region of interest 
(around 2997 keV [0]), an a-source was permanently placed close to the ZnSe crystal. The source 
consisted in an Uranium solution, covered with a thin mylar foil to absorb part of the as energy 
and to produce a continuum spectrum in the range 1-4 MeV. Finally, during a calibration run, an 
AmBe neutron source was placed close to the detector in order to produce high energy ys. 

The LD is calibrated with a 55 Fe source permanently faced to the LD surface opposite to the 
ZnSe. The source emits two X-rays at 5.9 and 6.5 keV. 

The main features of the detectors are summarized in Table |T} The rise and decay times of the 
pulses are defined as the time difference between the 90% and the 10% of the leading edge, and the 
time difference between the 30% and 90% of the trailing edge, respectively. The intrinsic energy 
resolution of the detector (Gbaseline) is estimated from the fluctuations of the detector baseline after 
the optimum filter application. Because of the better energy resolution, in the following we restrict 
our analysis to ZnSe-R. 



Table 1. Parameters of the bolometers. Amplitude of the signal before amplification (As), intrinsic energy 
resolution after the application of the optimum filter ((Jbaseiine), rise (z>) and decay (t</) times of the pulses. 
ZnSe-R is the thermistor chosen for the data analysis. 





As 


^baseline 








[jUV/MeV] 


[keV RMS] 


[ms] 


[ms] 


ZnSe-L 


10 


5.2 


5.6 ±0.3 


30 ±4 


ZnSe-R 


50 


2.4 


5.7 ± 0.2 


23 ±2 


LD 


980 


0.069 


2.8±0.1 


9.4 ± 0.5 




Figure 2. Detected Light vs Energy measured by the ZnSe in a calibration run. The energy axis has been 
calibrated using the most intense y peaks. The ZnSe crystal is faced to a smeared a source (upper band) 
and to an AmBe neutron source (high energy ys in the lower band). The events characterized by low energy 
(below ~ 1 MeV) and large detected light are particles that interacted in both the ZnSe and LD. The arrow 
points to the 82 Se Q-value (2997 keV). 



3. Scintillation 



The measured light as a function of the energy deposited in the ZnSe bolometer is shown in Fig. ||. 
Two regions can be clearly distinguished: the upper band, populated by a events provided by the 
smeared a source (continuum below ~ 4.5 MeV ee ) and by internal and surface a contaminations, 
and the lower band, populated by ys produced by the neutron source. The Light Yield of /3 /y 
events (LYj3 /y), defined as the amount of detected light per particle energy, does not depend on the 
energy. Fitting the /3/y band with a first order polynomial function, we obtain LY^ = 6.416 ± 
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0.008 keV/MeV. The Light Yield of as (LY a ) is larger than LY^/y, unlike in other known scintil- 



lating crystals [|28|]. This behavior, already observed in other ZnSe crystals operated at cryogenic 
temperatures, is not yet fully understood. Nevertheless, as it will be shown in Sec. |], this does not 
limit the discrimination capability. 

The continuum produced by the smeared a source was fitted with a first order polynomial 
function, resulting in ]_Y s ™ eared = (29.70 ± 0.17) keV/MeV. The fit revealed the presence of an 
energy threshold for the scintillation, which is E™ r ^ erf = (230 ± 12) keVa. At larger energies, 
where only a peaks due to bulk contaminations are present, we obtained LY„" //: = (26.62 ± 0.86) 
keV/MeV and E^ sh = (180 ± 140) keVa. The large error on the threshold is due to the large 
distance of the peaks from the origin, and does not allow a trustworthy comparison with the smeared 
source. However it is clear that LY b ^' Ik is lower than YX s ™ eared _ 

Understanding this discrepancy is not trivial. Indeed, for a given light detector and experimen- 
tal set-up, the light yield includes not only the scintillation process that produces light, but also the 
light transport from the luminescent centre to the detector. Even if we are not able to disentangle 
these two processes, we performed some tests to investigate the effects on the discrepancy between 
UYu' lk and YX s ™ eared induced by the crystal self-absorption, by the reflection on the mylar foil and 
by the energy dependence. In detail: 

• we checked whether the lower LY of bulk events could be due to the ZnSe self-absorption. 
This test was performed by placing the smeared source and the LD on opposite sides of the 
ZnSe, in such a way that the light emitted by as impinging on the surface had to travel across 
the entire crystal before reaching the LD. The results of this test are the ones already reported 
above and show that LY" ed > LY* ulk even if we expect the maximum re-absorption of the 
scintillation light emitted by the smeared as. 

• we investigated a possible energy dependence of the LY by facing a 224 Ra source emitting 
high energy as to the lateral surface of the ZnSe crystal, so that we could reproduce the 
same a particles from the bulk contamination and compare external and internal sources at 
the same energy. We obtained hYf ernal = (29.6 ± 0.1) keV/MeV from 4 to 6 MeV. This 
value is compatible with ]_Y™ eared pointing to the hypothesis that the larger LY is not due to 
energy-dependent effects. We underline that this test was performed without reflecting foil, 
to make sure that the larger LY of external as could not be attributed to the reflection on the 
mylar sheet. 

The experimental data show that in every test the light yield from bulk interactions is lower than the 
one produced by external events. Further measurements are needed to deepen the understanding 
of this behavior, that could be attributed to an effectively higher light production (the bulk of a real 
crystal is different in many aspects from the thin layer close to the surface), or to the different light 
collection for surface/bulk events. Indeed, it is well known that the light collection is particularly 



non uniform when dealing with cylindrical crystals [|29|]. 



To evaluate the discrimination capability at the OvDBD energy, we performed Gaussian fits to 
the light emitted in a and j3/y interactions, excluding the outliers due to a interactions in which a 
leakage of light is detected. From the fits we derived the mean value (/i) and the standard deviation 
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(a) of the light for both the a and j3/ 7 interactions. Since the discrimination capability increases 
with energy (as one can see in Fig. fy, we calculated /J. and a in several energy intervals and fitted 
the energy dependence of ^t(E) and <7(E) with polynomial functions. 
We defined the Discrimination Potential as a function of the energy as: 

\Ua(E) - llaJE)\ 

DP{E) = ^= ' mV 71 (3.1) 
ol{E) + oj y {E) 

and found DP = 17 at 2997 keV. It has to be remarked that this is just and indication of the 
capability of rejecting the a background. As it can be noticed from Fig. ||, there is a considerable 
number of as in which a large amount of light is lost. This behavior, likely due to surface effects, 
is particularly evident looking at the peak of 210 Po (~ 6.5 MeV ee ), a common contaminant found 
on the ZnSe surface. The loss of light from a particles constitutes a non-negligible and hard to 
estimate background to j8 //events. As it will be shown in Sec. the pulse shape of the light signal 
carries information on the particle type, irrespective of the amount of light collected, allowing a 
safe identification of all a events. 

We finally studied the LY at very low energies, which is of particular interest for experiments 
aiming at the detection of Dark Matter interactions. These experiments need to disentangle the 
signal produced by nuclear recoils below 30 keV (following the hypothesis that Dark Matter is 
made of WIMPs [ j30| , [3~l|]) from the background induced by j8/ys. We analyzed the recoils following 
the a decay of 210 Po. Since this isotope is deposited on the surface, the a particle can escape 
without releasing energy, while the nuclear recoil is absorbed in the crystal. We observe events 
centered at 139.6 ±0.6 keV ee in the ZnSe, i.e. 35% more energy than the nominal value (103 keV). 
The accuracy of the calibration function has been checked down to 511 keV, where it shows a 
deviation less than 1%. The extrapolation at lower energies is expected to maintain or reduce this 
deviation. From the fit in Fig. ||[ we evaluate the light emitted by nuclear recoils as < 14 eV at 90% 
C.L., corresponding to LY m . < 0.140 keV/MeV at 90% C.L. We evaluate the light emitted in the 
range 10 - 30 keV nr as < 1 - 4 eV at 90% C.L. for nuclear recoils, and 90 - 260 eV for j8 /7s. 

In this test the energy threshold was not optimized and set at 70 keV ee (~ 50 keV nr , assuming 
the 35% miscalibration), while to be competitive with present experiments it should be below 
10 keV nr . Given the baseline fluctuation (2.4 keV ee = 1.7 keV nr ) the required threshold could be 



reached with a trigger based on the optimum filter [32fl. However, to obtain a DP between j3/ys and 
nuclear recoils at least larger than 3, light detectors with baseline noise less than 20 eV RMS are 
needed. The light detectors currently being used are far from this value (see Tab. [I]) and obtaining 
resolutions better than 70 eV RMS does not seem to be achievable with the present technology. To 
search for Dark Matter interactions in ZnSe bolometers, new light detection technologies must be 
introduced. 



4. Energy Resolution 

The energy resolution on the OvDBD signal is estimated from 7 lines, which produce the same 
bolometric and scintillation response of j3 particles. The energy resolution is found to be 13.4 ± 
1.0 keV FWHM at 1461 keV ( 40 K contamination in the cryostat) and 16.3 ± 1.5 keV FWHM at 
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Energy [keVee] Detected Light [keV] 

Figure 3. Left: Energy region around the nuclear recoils from the 210 Po decay. The detected energy of the 
recoils, using the j3 //calibration function, is 35% larger than the nominal value (103 keV, see discussion in 
the text). Right: Distribution of the light emitted from the nuclear recoils, selecting the events marked in red 
in the left figure and applying the cut Detected Light < 0.3 keV to remove the j3 jy background. 
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Figure 4. Left: Correlation between the detected light and the energy released in the ZnSe for 2615 keV y 
events. Right: Energy distribution in the ZnSe before (gray) and after (black line) the combination with the 
light signal for the same events. 



2615 keV ( 2US T1 from calibration). These values are worse than the baseline resolution, which is 
5.6 keV FWHM (see Tab. [I]). Looking at the detected light, we observe that there is a positive 
correlation with the energy measured in the ZnSe (Fig. || left). This is in contrast with other 
scintillating crystals, where the correlation is negative as expected from the conservation of energy 
[]|]. Presently we are unable to explain this behavior, however we can take advantage from the 
correlation to improve the energy resolution. 

The light L and energy E in the ZnSe are two correlated estimates of the same observable, i.e. 
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the particle energy. Defining L 1 = L/LYp/y, we write the combined energy variable as: 

E comb = {l-w)E + wL' . (4.1) 
The weight w that minimizes the variance of E comb is: 



ol(l-po LI /o E ) 
ol + ol,-2po E o v 



where aj and o]j are the variances of E and L' , respectively, and p is the correlation between them. 
The expected variance of E comh is then: 

2 _ P 2 ) 

o E + o[,-2po E o L > 

The application of this algorithm to data improves the energy resolution considerably, in partic- 
ular at high energies where the energy resolution at 2615 keV improves from 16 to 13 keV FWHM 
(see Tab. || and Fig. fright). 

The resolutions of the light detector and of the ZnSe, as well as their correlation, depend on the 



energy. To apply Eq. |4JJ to the entire energy spectrum, even where there are no peaks, a specific 
weight should be used (Eq. fO| ) at each energy. Since the resolutions and the correlation are slow 
functions of the energy, w(E) has been estimated from polynomial functions fitting the energy 
dependence of a E , o E > and p. 

Table 2. Energy resolution before and after the combination with the detected light. 





ZnSe 


ZnSe and Light 




[keV FWHM] 


[keV FWHM] 


1461 keV 


13.4 ± 1.0 


12.2 ±0.8 


2615 keV 


16.3 ± 1.5 


13.4 ± 1.3 



5. Pulse shape discrimination 



Given the different light and heat yields of a and /3/y interactions, we tested whether also the 
shape of the signals carries information on the type of interacting particle [33, [34] ]. 

The average pulse of 2615 keV y events in the ZnSe and the corresponding average scintil- 
lation pulse are shown in Fig. ||. The pulses are fitted with a model developed for Te02 bolome- 
ters [ p5| , [36[ ]. The model originally included the thermistor and electronics responses and the re- 
sponse of the bolometer, which describes the time development of the thermal phonon signal after 
an instantaneous energy absorption: 



AT(0 



. e -t/* + ae-'l^ + (1 - a)e-' /r " 2 



(5.1) 



where A is the thermal amplitude, x r the rise time, two decay constants and a is a weight 
(0 < a < 1). In order to introduce the model in this work, the above equation has been modified to 
include energy releases that can not be considered instantaneous, such as the scintillation process. 
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Figure 5. Average pulse of 2615 keV y events in the ZnSe (left) and the corresponding average scintillation 
pulse (right). The dashed lines represent the fast and slow components of the fit function. 



A function which has been found to well reproduce the observed signals, both in the ZnSe and 
in the light detector, is: 



AT(t) 



- e -t/* + ae-tfa + (1 _ a )e-'l^ 



+A 2 



- e -t/*r + e -t/%B 



(5.2) 



In Fig. |5| the total fit function and the two components are displayed. In the ZnSe the amplitude of 
the fast component is given by A\, while the slow one by A%. In the light detector the amplitude of 
the fast component is A%, while the slow one is A\. We suppose that the fast decay of the signal can 
be attributed to phonons, while the slow one is affected by the scintillation mechanism, which is 
different for a and j3/y interactions [|37|]. Applying the fit to each event, we noticed that the ratio 
of the slow and fast amplitude components (A s /Af) is related to the particle type. However the fit 
is too sensitive to the detector noise and to minimization problems and cannot be used to estimate 
a reliable discrimination parameter. 

Following Refs. [38, |39|], we developed a bi-component optimum filter algorithm to lower 
the noise and improve the reliability of the parameter A s /Af. To explain the algorithm, we first 
remind to the reader the working principles of the standard optimum filter. Given a signal + noise 
waveform f{t) = A- S(t — to) +n(t), the best estimate of the amplitude A, under the hypothesis that 
the noise is stationary, is obtained from the x 2 of the residuals in the frequency domain: 



\f(co)-A-S(co)e- 
N(co) 



imt 



(5.3) 



where S(t) is the ideal pulse, estimated as the average pulse, and N(co) is the noise power spectrum 
of the detector. The parameter to accounts for any possible jitter between the observed signal and 



S(t). Minimizing Eq. 5.3 with respect to A, we obtain: 



A{t Q ) = hY J 



S*((0)e 
N((0) 



l(Ot 



-m 



(5.4) 



where h = [£ ffl \S((o)\ 2 /N((o)] l . To estimate A, Eq. ||3] has to be minimized also with respect to 
to. This however cannot be done in an analytical form. It can be demonstrated that the minimum of 



Eq. 53 with respect to to and A is equivalent to the maximum of A(to) . Since Eq. 5^ is an inverse 
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Fourier transform on to, A(to) can be seen as the filtered signal in the time domain. In this domain, 
the best estimate of A (A) can be extracted with a maximum search algorithm. 

In the bi-component filter the waveform is decomposed as f(t) = A s S s (t — to) +AfSf(t — to) + 
n(t), where in our case S s (t) and S/(t) are the ideal slow and fast components of the signal. The 
minimization of the % 2 with respect to A s and Ay gives: 

Ut )\ = (l\s s \ 2 n- 1 Us:s f )N- l Y l (zsy^N-'A 
Mto)J \Us* s s f )N- 1 L\s f \ 2 N- 1 ) ' [zsy^N-'f) 

where the sums run over CO and S s ,Sf,N and / are to be intended in the frequency domain. In this 
case the minimum of the % 2 with respect to to does not correspond to the maxima of A s (to) and 
Ay (?o)- Therefore we estimate to by scanning the % 2 around to = and then evaluate Ay = Ay(?o) 
and A* =A s (t ). 

The application of the bi-component optimum filter is shown in Fig. ^ using the fit components 
in Fig. H as estimates of S s j(t). As the figure shows, the separation is very evident in the light 
detector and less evident in the ZnSe. In the case of the ZnSe we estimate DP = 2, while for the 
light detector DP = 11. As it will be shown in the next section, the pulse shape selection applied 
to the light detector allows one to tag a events with light yield compatible to that of j3/ ys, thus 
dramatically reducing the background. 



6. Background 

We performed a background run of 524 hours in order to evaluate the internal contaminations of 
the crystal. To further decrease the environmental j8 /y background, the crystal was shielded with 
ancient Roman lead, featuring an activity lower than 4 mBq/kg in 210 Pb j40|] - Figure [7] shows the 
detected light as a function of the particle energy measured by the ZnSe. We notice that the pulse 
shape cuts allow to completely identify the interacting particle. 

In Figure ^ (left) the distribution of j8/y events is reported. We can identify the decay of 75 Se 
(Tj/2 = 119.779 days, Q = 863.6 keV). 75 Se is produced via 74 Se neutron capture and decays via 
electron-capture (100%) with a complex combination of de-excitation ys and X-rays, producing 
a peak at ~ 410 keV and a series of peaks between 150 and 300 keV (Fig. ||). We analyzed the 
decay rate of the 410 keV peak as a function of time, obtaining T!/ 2 ( 75 Se) = 102 ± 18 days, which 
confirms our expectations. Another isotope that is produced by neutron activation in ZnSe is 65 Zn 
(Tj/2 = 244 days, Q = 1359.9 keV). This isotope decays via electron capture and produces a peak 
at about 1350 keV. These contaminants do not contribute to the OvDBD background of 82 Se, 
because of the low Q- value and the short half-life. The other events in the j8 /y spectrum can be 
attributed to 40 K and 232 Th contaminations of the environment. 

Finally, in the inset of Fig. || (left), the zoom on the energy region of interest of 82 Se is reported. 
One can see that only one event above 2615 keV is observed. This event, however, is in coincidence 
with high energy y's seen by other scintillating crystals that we were testing in the same set-up. We 
believe that it can be ascribed to a muon interaction in the surrounding materials. This kind of 
events can be easily suppressed in bolometric arrays by applying time-coincidence cuts, as done in 
CUORICINO ^ or, ultimately, by surrounding the cryogenic facility with a muon veto. 
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Figure 6. Left: Ratio of the amplitudes estimated by the bi-component optimum filter algorithm as a function 
of the energy released in the ZnSe. The events tagged as a are marked with red, while j3 /7s with blue. 
Right: Histogram of the ratio. Heat in ZnSe in the top row, light in the LD in the bottom row. The a and 
j3 /y populations are well identified in the light detector but less in the ZnSe. 



In Fig. ^ (right) the a background is reported after cutting on both pulse shape and energy of 
the light channel. The activities of the most intense peaks, reported in Table [3j show that 232 Th and 
238 U are in equilibrium with their daughters. 

In principle, a large internal contamination in U could be worrisome because of one of its 
daughters, 214 Bi, which /3 -decays with a Q- value of 3272 keV. However, 214 Bi decays with a BR 
of 99.98% in 214 Po, which a-decays with a Q-value of 7.8 MeV. The half-life of 214 Po, which is 



Table 3. Activity of the isotopes belonging to Th and 2J8 U chains. 



Chain 


Nuclide 


Activity 
[MBq/kg] 


232 Th 


232 Th 


17.2 ± 4.6 




228 Th 


11.1 ± 3.7 


238TJ 


238tj 


24.6 ± 5.5 




234 u 


17.8 ± 3.3 




23() Th 


24.6 ± 5.5 




226 Ra 


17.8 ± 3.3 




210p o 


< 90.9 ±10.6 
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Figure 7. Data from 524 hours of background runs. The different colors mark j3/y events (blue) and a 
events (red) selected with pulse shape cuts on the light channel. The arrow points to the Q-value of the 82 Se 
decay. The gray dots indicate particles interacting in both the ZnSe and the light detector (double hit events), 
and events with no detected light (dark events). Double hit events are identified through their pulse shape in 
the LD, while dark ones through the shape in the ZnSe. 
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Figure 8. Simulation of the 75 Se decay in the ZnSe. Energy resolution set to 7 keV FWHM. 
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Figure 9. 524 h background run. Spectrum of /3 /y (left) and a events (right). Data are calibrated using y 
sources, as result miscalibrated by +22% with respect to their nominal energy. 



about 160 jits, is extremely short compared to the time development of bolometric signals, which 
is hundreds of ms. For this reason, the /3/y emission of 214 Bi is simultaneously followed by a 7.8 
MeV a particle and does not represent a problem for the background, as the superimposition of the 
two signals lies at much higher energies than the region of interest. 

When evaluating the background due to 214 Bi we considered also its decay in 210 T1, which 
occurs via a emission with a BR of 0.02 %. This contribution, however, can be easily suppressed: 
210 T1 /3-decays with a Q-value of 5489 keV and a half-life of 1.30 min. The background induced 
by the j3 emission can be rejected using the delayed coincidences with the a of 214 Bi. In this case, 
the dead time induced by the delayed coincidences is almost negligible, as the half-life of 210 T1 
is very short and the BR for this decay is extremely small. As a consequence, even an internal 
contamination in U of 25 /iBq/kg would produce a background in the energy region of interest 
of only 3 x 10 -4 counts/keV/kg/y. 

A more dangerous background source is the one due to the internal contamination of 208 T1, 
that belongs to the 232 Th chain and /3-decays with a Q-value of 5001 keV. 

Considering a 232 Th contamination of 17 /iBq, we expect a background in the energy region 
of interest due to 208 T1 of about 2.5 x 10 3 counts/keV/kg/y, which represents the ultimate limit for 
the background. The background due to 208 T1 could be suppressed using the coincidences with 
its parent, 212 Bi, that a-decays with a half-life of 3 minutes. In this case, however, the introduced 
dead-time must be considered. A background reduction factor of 3, for example, requires a delayed 
coincidence of 10 min, which would imply a dead time of about 10%. 

Summarizing, the internal contaminations of 214 Bi and 208 T1 do not represent a dangerous 
problem for the achievement of a background of 10~ 3 counts/keV/kg/y, which is the goal of the 
LUCIFER experiment. In any case, purification techniques are under study in order to reduce bulk 
contaminations of the crystals and, therefore, the dead time induced by delayed coincidences. 
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7. Conclusions 

In this paper we analyzed the performances of a 431 g ZnSe crystal operated as scintillating 
bolometer. We measured the energy resolution at 2615 keV, close to the OvDBD of 82 Se (2997 keV). 
By combining the light and heat signals, the resolution has been estimated as 13 keV FWHM. 

We presented the results of a 524 hours background run to identify the internal radioactive 
contaminations of the crystal. We measured 232 Th and U contaminations of the order of tens of 
/iBq/kg, which is compatible with the low background requirements of a OvDBD experiment. In 
this run, no event passed the data analysis cuts against the background in the energy region of the 
decay, showing the potential of this detection technique. 

We analyzed the scintillation properties of the detector and we demonstrated that the infor- 
mation carried by the light signal allows to identify the nature of the interacting particles in two 
different ways. Exploiting the different light yield of as and j6/ys, one can reject the background 
due to a particles within the energy region of interest for the 82 Se decay. We developed an algo- 
rithm, based on the shape of the light signals, to tag a particles with light yield compatible with 
j3/y ones. 

Finally, the good identification of nuclear recoils at 100 keV allowed us to set the requirements 
for Dark Matter searches with ZnSe bolometric detectors. We estimated that nuclear recoils of 
210 Po are detected with an amount of energy 35% greater than the nominal value, using a calibration 
function estimated from y lines. This value has to be confirmed with dedicated tests, together 
with the evaluation of the calibration function of Zn and Se recoils. By using light detectors with 
baseline noise lower than 20 eV RMS, we estimate that j3/y events could be discriminated from 
nuclear recoils events of 10 keV, thus setting indicative requirements for Dark Matter searches in 
LUCIFER. 
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